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The solution of many problems of statistical data processing is associated with the generalization and 
identification of patterns in the processes under study based on observations of their real flows or physical 
experiments. These processes are always influenced by many random and discrete factors. In a physical 
experiment, statistical data are collected, which are usually incomplete. In such cases, the only research 
method is the mathematical analysis of the information received, where heuristic self-organization methods 
are widely used [1]. One of them is the method of group accounting for an argument (GAM), which is based 
on a multi-row selection of particular models with successive complications. 


A.G. Ivakhnenko [1] gives a list of programs based on GAM algorithms. The logical substantiation and a 
fairly accurate presentation of the principle of self-organization contributed to the deepening of these 
studies. 


This article is devoted to the development of algorithms for parametric identification of the qualitative and 
dynamic states of the system, respectively. 


Determination of qualitative estimates of geophysical parameters under various artificial impacts on the 
reservoir system is one of the most important issues in the analysis of oil and gas field development 
processes. To identify the information content of features, a generalized parametric identification algorithm 
has been developed to study the dynamic states of the system based on the principle of “Heuristic self- 
organization”. 


The input data of the object of study here are the results of a physical experiment, which are accepted as 
training information, and the information obtained with the help of the constructed model is as a verification 
one. 
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To select the optimal object model, the algorithm provides for checking the standard deviations using the 
“left corner” principle. 


Main part 


Let us consider the construction of heuristic models of multipara meter objects based on the results of a few 
physical experiments. Let some information about the studied process be collected from observations with 
volume @, which is the starting point for building the model. This volume of information is divided into 
three parts: training @), verification @, and examining @,. In contrast to the algorithms of the method of 
group accounting for an argument (GAM), this separation is not carried out in a unique way. The model built 
by the GMDH method is optimized by the ranking of the initial information @. So there is a third possibility 
of choosing the optimal model. For this, root mean square errors are used. @,, calculated by the optimal 
model of each ranking of the initial information@. 

Thus, the task of constructing a mathematical model of multi-parameter processes based on the MGAA 
reduces to the following: incomplete information is given @, consisting of system parameter values 
I E Xij > 
x,;. Required to define a dependency function y = f (xX, X,,...,X,) and its predicted value y, (n -number 


i=1,n, corresponding to the studied parameter y j» J=1,R and parameters at the predictive point 


of measured parameters in the system, R - number of observations). 
Stages of building a three-criteria selection model. 


1. Transition from dimensional to dimensionless variable formulas 
(1) 


X. 


‘max? Ymax - Maximum values. This stage makes it possible to calculate the values necessary for constructing 
the model, regardless of the setting of the system of parameter dimensions, therefore, all parameters can be 
represented in the dimensions accepted within the studied limits. This makes it easier to prepare the initial 


information. 


2. Pre-processing of initial information. We calculate the range of parameters in the response function, the 
mean value of the variance and the standard deviation using the formulas. 


F = Ximax ~ Ximin 

Vie Vine 

a R 

x= Y x, /R, (2) 


jal 
yO a 
S? =a, - x) /(R-D, S =S, = LN. 
j=l 


3. Selection of the examining sequence from the initial information. To do this, we calculate the distance 
between the point and the steel points of the definition using the formula 


[Nn — 
Pi = Yj — Hu) F=LR (3) 
i=l 


Source Information points @ sorted by increasing distance p}, < P, <...< Op. Points corresponding to the 


first M values 9;, form examining sequences @,, the rest are under study a, . 
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(ay) 


4. Separation Ou on the ® and © 


with volume R, and R, For this point from Oy , located in 
ascending distance from the forecast point, are collected in various combinations and are formed > and 
On For example, ®. consists of points with even numbers, a @n_ odd. Vice versa, @n consist of even, 


@._ from the odd ones. The first two in 2 , the third is in On etc. Can arrange M, combinations. 


5. Pairing Ak and X, parameters and construction of particular descriptions. Dependencies can be taken 
from support functions 


Y, = Ay + A,X, + A,X, + 4X, X; (4) 
or 
_ 2 2 
Y, = Ay +A,X, + A,X, + A,X,X, + A,X, + AsX, 
kK=1,N-1, 1=2,N, 1=2,M,, 
where M, =1;, - number of private descriptions. 


Dependencies of the form (4) are constructed by the least squares method in @, . The first selection criterion 
is used 


S £52 . 
AN) + min. 


1 j=l 
Here Yy calculation according to dependence (4) of the value of the studied parameter y ;. 


6. Calculation of quantities 


R 


1 *\2 
Ô; zO 


2 J= 
and ordering particular models in ascending order 6? 


ô < 07 <... < 8H, 


After that, from the set of particular descriptions, the first. 
Let's move on to the next set of selections. The selection process continues until 


+3 
pays 


«2 
0 >é or 01 >e db; > 6, (5) 


bo 


Here č required calculation accuracy; 6,” mean square error values from the previous selection 
corresponding to the best of particular descriptions. 
+2 

Condition (5) corresponds to the “left corner” principle. After fulfilling the condition Oo, > 61. The next 
selection does not improve the model; in the general case, condition (5) for this statistic is also satisfied for 
all subsequent selections, i.e. further improvement of the model is not possible. 

This is how the only optimal model is determined that satisfies the criterion of minimum mean square error 
and corresponding to M - th ranking of the initial information ©, 


uc * 


7. Calculation of the value of the root-mean-square error for the examining sequence in the optimal model 
using the formula 
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M, 2 


a yyy). 
a my y’) 


where y? - calculated value corresponding to the table y ; and calculated by the optimal model. 


8. Repeat steps 5-7 for another ranking. This determines the set of the best models corresponding to the 
ranking options. 


9. Model corresponding to the minimum of the set 8? , is considered the best and its accuracy is determined 
by the formula 
62, = minô? 
Thus, the best model of multiparametric processes based on the principle of self-organization is built. It is 
practically impossible to get the model form, as in [1], because the addition of ranking options will again 
increase the number of terms in the polynomial. Apparently, this is not necessary. If you want to get the 
numerical value of the forecast, this can be achieved by calling the function procedures. 


When using the above algorithm, its features should be taken into account. If one of the pairs of arguments 
X,, X, takes uniform values, i.e. they can be represented as 


Xi, = jh, x), = jh,, j=LR, 


then the determinant of the normal system of equations is equal to zero. In this case, it is impossible to 
obtain solutions to the normal system of equations and, consequently, to build a model for this pair. We 
exclude the model corresponding to this case from the set of models and proceed to the next pair of 
parameters. Therefore, the total number of models in each row of selections may be less than C Fe : 


Consider the problem of determining the parameters of the system if the value of the response function is 
given. Let the best model of the system be built, it is required to determine the values of the parameters 


(x), X9,..., X?) under which the function f (x), x},..., X?) = fy. It is known that this problem reduces to 
finding the minimum of the functional 


F(x, X NR et ss BD Yo 
with restrictions on parameters X min SX; < Xpmax - 


An algorithm and a program for determining the global minimum of the functional based on the random 
search method were developed [2]. The program can be successfully implemented in this case as well, for 
which the GAM function procedure should be used as a function. 


The proposed MGAA algorithm for constructing a model of multi-parameter processes is implemented on a 
computer in the C++ language. We use the proposed technique to simulate the well killing process, 
considering it in two stages. 


1. The process of filtering the washing liquid. The permeability of the reservoir depends on the 
physicochemical properties of the drilling fluid and the properties of the reservoir. The following 
dependency is used: 


K im = Kim, V, u,y,B, C K,,), (6) 


c? 


where K, - initial air permeability of the core to saturation with solution; K „- permeability of the core with 
washed crust; T „- duration of core saturation with solution; V - the volume of the passed filtrate; 4- 
solution viscosity;  - specific gravity of the solution; B - water yield in 30 min. solution; C,- static shear 
stress. 
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The initial information is given in table.1. 


Table 1 
K, To V H 7 B C K,,, 
761,4 18 0,1 117 1,02 7 18 168,3 
698,2 16 8,6 63 1,02 38 0 150,4 
1037,3 19,5 4,9 36 1,04 8 84 507,8 
1025,9 20 9,4 50 1,05 8 12 428,2 
593,5 14 8 28 1,062 9 0 128,2 
698,2 14,5 3 28 1,07 8 0 287,0 
1047,3 16 3,9 56 1,12 10 24 204,9 
1044,5 17,5 11 38 1,13 18 20 448,1 
698,2 17 9 28 1,14 12 9 229,0 
530,1 20 4,5 28 1,13 32,5 0 198,3 
Table 2 
K,,, To H } B C. K, 

168.3 2,5 117 1,02 7 18 725,5 

150,4 3 62 1,02 38 0 553,9 

507,8 1,75 36 1,04 8 84 946,01 

428,2 1,66 50 1,05 8 42 814,56 

128,2 3,5 28 1,062 9 0 539,5 

287 2,75 28 1,07 8 0 584,4 

204,9 4 56 1,12 10 24 827,37 

448,1 2,33 38 1,13 18 20 857,33 

229 2,66 28 1,14 12 9 609,53 

198,3 2,75 28 1,13 32,5 0 416,66 


2. The process of displacement of the filtrate by gas. In this case, the permeability of the reservoir can be 
slightly improved. The dependence is considered 


K,=K (Kk Is H, y» B, Cae) 


HM? 


(7) 


where K,- core permeability after blowing; T 


y 


„p` air purge time. 


Initial information is given in tables 1, 2. 


b 


HM 


The proposed method was used to construct dependences (6) and (7). Where K’ and K? - calculated 
values of the parameters according to the optimal model. 


The change in the initial permeability of the rock is judged by the following data (Table 3): 
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Table 3 
Kin K”, (K nu _ K?) / K?, 
168,3 176,96 -0,05140 
150,4 162,48 -0,0732 
507,8 502,67 0,0101 
428,2 437,26 -0,0211 
128,2 133,88 -0,0443 
287,0 264,51 0,0783 
204,9 201,31 0,0175 
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448,1 458,11 -0,0223 
229,0 249,29 -0,0886 
198,3 191,51 0,0342 


The change in rock permeability after gas blowing is characterized by the following data (Table 4): 


Table 4 
K, K? (K, -KIK 
725,5 703,86 0,0298 
553,8 548,69 0,0092 
946,01 931,18 0,01577 
814,56 827,46 -0,0158 
539,5 506,39 0,0614 
584,4 567,23 0,0294 
827,37 828,03 -0,0008 
857,53 874,02 -0,0198 
609,53 634,77 -0,0411 
416,66 445,15 -0,0684 


The optimal model was obtained for the first stage after 27, for the second after 6 rows of selections. The 
mean square errors are 0.0796 and 0.00034. The reason for stopping the selection of models was the “left 
corner” principle for both stages, i.e. improvement of the mean square error in the future is impossible. The 
relative error ranges from 1 to 8%. Similar results were obtained with other data of the problem. 


Naturally, in the cases described, it is fundamentally impossible to bring the model into an analytical 
connection of the response function of the studied parameters. Such a model would include at least 28 terms 
of varying degrees of variables. 


Thus, the proposed method can be used to build an optimal model for changing the filtration properties of 
reservoirs when flushing fluids with different properties move through them [4]. 


A good agreement between the experimental and calculated data shows the acceptability and effectiveness of 
the proposed parametric identification algorithm when developing a formulation of lightweight drilling 
fluids with and without solids to improve the technological and technical and economic development 
indicators when drilling wells in difficult conditions, as well as when killing gas wells on the period of 
production of repair work at the fields 
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